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For a quantum-mechanically spread-out particle we investigate a method for determining its 
arrival time at a specific location. The procedure is based on the emission of a first photon from 
a two-level system moving into a laser-illuminated region. The resulting temporal distribution 
is explicitly calculated for the one-dimensional case and compared with axiomatically proposed 
expressions. As a main result we show that by means of a deconvolution one obtains the well known 
quantum mechanical probability flux of the particle at the location as a limiting distribution. 
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CN . I. INTRODUCTION 

■ An important open problem in quantum theory is the question of how to formulate the notion of "arrival time" 
of a particle, such as an atom, at a given location, i.e. the time instant of its first detection there. This is clearly a 
very physical question, but when the extension and spreading of the wave packet is taken into account, a satisfactory 
formulation is far from obvious. The problem of time in quantum mechanics, both for time instants and time durations 

, such as dwell time, has received a great deal of theoretical attention recently |^ . When the translational motion of 

' the particle can be treated classically, a full quantum analysis of arrival time is in fact not necessary. This is the case for 

[ fast particles, and therefore arrival times are presently measured mostly by means of time-of-flight techniques, whose 

■ analysis is carried out in terms of classical mechanics. Problems, though, arise for slow particles for which the finite 
^\ I extent of the wavefunction and its spreading become relevant, such as for cooled atoms dropping out of a trap. Then 

. a quantum theoretical point of view is needed. It is therefore important to find out when the classical approximations 
' fail and to propose measurement procedures for arrival times in the quantum case. Since the theoretical definition of a 
: quantum arrival time is still subject to debate it is necessary to determine what exactly such measurement procedures 
' are measuring and to compare such operational approaches with the existing, more abstract and axiomatic, theories. 
Oh! -A- simple one-dimensional example is the arrival time at a; = of a particle described by a wave packet moving 
' to the right. The probability of finding it still on the left side at time t is J^^dx\'ip{x,t)\'^ . Then it would seem 
^ [ natural to assume that the probability of arrival in dt is given by the decrease in dt of this integral, i.e. that the 
^ • arrival-time probability density is its negative derivative. Since p + dJ^/dx = 0, where is the usual quantum 
mechanical probability current, one immediately finds that the arrival-time probability density should be given by 
• • I J^{0,t). For an ensemble of classical particles the arrival time distribution is also given by the flux, so everything 
seems to fit nicely. However, the quantum flux for a wavefunction formed entirely of positive-momentum components 
may become negative at certain times j^. As a way out Leavens Q has proposed, following a Bohm trajectory 
^ ■ analysis, to use its normalized absolute value. 

The difficulties to formulate a quantum arrival time concept were posed most prominently by AUcock j5| , and there 
are several attempts to overcome these [Q, ||. Kijowski [Q, in particular, obtained a time-of-arrival distribution for 
free motion from a set of axioms modeled after the classical case; for a general investigation see Werner The 
distribution of Ref. ||^ has been studied, compared to other approaches, and generalized by some of us for systems 
subject to interaction potentials and to multi-particle systems j9[ ^ ^ ^ ^ 

No procedure how to measure the proposed distributions was given, nor is one known today. The gap between 
experiment and axiomatically defined quantities has been commented on and considered to be worrying by Wigner 
and others; for a review see Ref. §. In this vein, several "toy models" for arrival-time measurements have been 
put forward by Aharonov, Oppenheim, Reznik, Popescu and Unruh (cf. also ^ ^^), but these models do not 
incorporate the basic irreversibility inherent in any measurement process. Irreversibility has been included by Halliwell 
[0 in a model based on a two-level detector in which the initial excited level decays due to the presence of the particle. 
However, the model remains somewhat abstract since no connection is made with any specific measuring system. 

An experimentally very natural approach to determine the arrival time of an atom is by quantum optical means. 
A region of space may be illuminated by a laser and upon entering the region an atom will start emitting photons. 
The first photon emission can be taken as a measure of the arrival time of the atom in that region. This approach has 
been proposed by three of us and Baute |l8[ , and in a preliminary study of the one-dimensional case a surprisingly 
good numerical agreement with Kijowski's axiomatic distribution was found in some special examples. 

From a fundamental point of view, however, immediate objections may be raised to this experimental, or operational. 
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approach. First, depending on the decay rate of the excited atomic level, the photon emission will not be instantaneous 
but will take some time, thus leading to a delay compared to some "ideal" arrival time of the atom. Second, the 
laser takes some time to pump the atom from its ground-state to an excited state, and therefore this also leads to a 
delay. Conceivably, the second objection might be overcome by progressively increasing the laser intensity, and the 
first objection by considering shorter life times so that in a theoretical limit one would arrive at an "ideal" quantum 
arrival time without the above shortcomings. Attractive though this seems at first sight, it does not work, as will be 
shown in this paper. The reason is a further difficulty - reflection. Although the laser couples only to the internal 
degrees of the atom, it will be seen that there is a nonzero probability for the atom to be reflected from the laser 
region without ever emitting a photon. Nevertheless, there is a way out of these difficulties, with a surprising result. 
The idea is to "subtract" the delays from the first-photon probability density by means of a deconvolution with an 
atom at rest. This results in a distribution which, for shorter and shorter life time of the atomic level, converges to 
an unexpected distribution - namely to J^, the quantum mechanical probability flux. The probability distribution 
for the first photon is non-negative and the emergence of possible small negative values is due to the deconvolution 
procedure. This connection to opens a way, to our knowledge for the first time, to measure the quantum mechanical 
probability flux. 

For simplicity, this paper considers only the one-dimensional case. The probability density for the emission of the 
first photon from a moving atom is calculated explicitly by means of the quantum jump approach [p^ . It is shown 
that large laser intensities lead to a large reflection probability. This in turn leads to a large non-emission probability 
and a first-photon probability density not normalized to 1. Then the problem of reflection versus time delay is 
discussed. Reducing the laser intensity leads to a pumping delay. It is shown that trying to reduce the emission delay 
by shortening the level lifetime leads in the limit to a free wave packet in the ground-state with no emissions. The 
delays are then removed by a deconvolution, and we discuss for which parameters the resulting expression is close 
to its "ideal" limit J^. For more practical purposes it is also shown that for a certain domain of parameters, which 
include those used in Ref. , the non-deconvoluted first-photon probability density gives a good approximation to 
and to Kijowski's axiomatic arrival time distribution. However, it is also pointed out that Kijowski's distribution 
cannot be obtained in a simple direct way as an exact limit of our operational approach. 



II. THE PROBABILITY DENSITY FOR THE FIRST PHOTON 



The Hamiltonian of a two-level atom of mass m, interacting with the quantized electromagnetic field E and a laser 
with (classical) field El is, in the usual dipole approximation and in the Schrodinger picture, 

H^-^+HA + Hp + eT>-{^ + ^L) (1) 
2m 

where 

Ha = i?.c.{|2)(2|-|l)(l|} (2) 



Hp = ^tC^kflkA^kA 



kA 

D = di2|l)(2| +h.c. 

with di2 the transition dipole moment between the states |1) and |2). It is assumed that the laser illuminates a half 
space, xi > 0, say. 

Let an atomic state |\E'(to)) be prepared at time ^o- By means of the quantum jump approach the atomic time 
development until the first photon detection is given by a (non-hermitean) "conditional" Hamiltonian 

\^{t)) =e-*^^(*-*°)/'"'|*(to)), (3) 

with the photon part traced away. In the interaction picture with respect to the internal Hamiltonian Ha one has in 
the usual rotating wave approximation 

H, = pV2m + e(xi) {|2) (I|e^'^-- + h.c.} - '-h^\2){2\ , (4) 

where the Rabi frequency cx d • E^^-* plays the role of a laser-atom coupling constant (with E^"-* the laser amplitude) , 
ki is the laser wave vector, and where 7 is the Einstein coefficient of level 2, i.e. its decay rate or inverse life time. 
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One can show that Eq. includes the Doppler effect, i.e. the laser driving depends on the atomic velocity through 
a frequency shift. The probability, N^, of no photon detection from to up to time t is given by [^9| 

Nt^\\e-^"^(*-'°^/^mo)W, (5) 
and the probability density, n(t), for the first photon detection by 

dNt , , 

n(t)^-^. (6) 

For simplicity, we only consider the corresponding one-dimensional problem, with the laser perpendicular to the 
atomic motion, so that the Doppler effect plays no role. With |1) = (p) and |2) = ("j*) the conditional Hamiltonian 
becomes, in matrix form 

To obtain the time development of a general wave packet under we first solve the eigenvalue equation 

{x) 



Hc^ = E^, where = ( ^,3)!^^ )• (8) 



Since for a; < there is no laser one obtains for 0^*^ the equations 

2m 2m 

Hence 0^^-' is a combination of e**^^ and e"**^'^ with k satisfying 

E = h^k^/2m = Ek (10) 

while 0*^^) is a combination of e*''^ and e~**^ with q satisfying 

E + ih-//2^h^q^/2m . (11) 

We now look for eigenstates of which correspond to a ground-state plane wave coming in from the left. Then k as 
well as E must be real, by boundedness, and for x < the eigenstate is of the form 



1 



^ikx _|_ p ^ — ikx 



with Im q > for boundedness, while Ri and R2 are refiection amplitudes yet to be determined. Note that although 
E = Ek is real the complete wave functions will not be orthogonal, in accordance with the non-hermiticity of He- 

To obtain the form of ^ki^) for x > we denote by |A-|-) and |A_) the eigenstates of the matrix ^ (^^ 

corresponding to the eigenvalues A±. One easily finds 



A± = --7±-V7^-4f]2 (13) 
|A±> = • (14) 



\2\±/n 

We exclude, for the moment, the limiting case A+ = A_. Note that |A±) are not orthogonal and have not been 
normalized. For x > 0, one can write as a superposition of |A±) in the form 



V2^*fc(x) = C+|A+)e''=+^ + C_|A_)e*'''-^, a; > 0. (15) 
From the eigenvalue equation iJc^fe — Ek^k, together with Ek = h^k"^ /2m, one obtains for a; > 



fc2 = fc2 _ 2m\±/n = k^ + imj/2h T im^/-f^ ~ Aft^ /2h, (16) 
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with lmk± > for boundedness. From the continuity of $fe(x) at a; = one gets, with |1) = (J) and |2) = (°), 

l + Ri = C+(1|A+) +C_(1|A_) 
i?2 = C+(2|A+>+C7_(2|A_). 

Similar relations result from the continuity of ^'i^{x) at x = 0, yielding 

C+ = -2fc(g + fc_)A_/L>, (17) 
C- = 2k{q + k+)\+/D, 
i?2 = fc(fc_ - k+)n/D, 

Rl = [X+{q + k+){k-k^)-\^{q + k^){k-k+)]/D, 
where the common denominator D is given by 

D = {k + k-){q + k+)\+ -{k + k+){q + fc_)A_ . (18) 
Thus Eq. ( p^ ) becomes, in components and for x > 0, 

</)«(x) ^ -^i^{(g + fc_)A_e''^-+--(g + fc+)A+e^'^--}, (19) 



^i^)(x) = ^2-{(g + fc_)e^'=^--(g + fc+)e^^-} 



The case A_ = A+ is obtained from this by taking limits. For later purposes we also consider increasingly large 7, the 
other parameters kept fixed. This leads to 

A+ w -^17727^0, (20) 

A_ w -i-f/2 + in^/2j-^ -ij/2, (21) 

kl ^ k^ + imfl^/h-f, k+ k, (22) 

fc^ « k'^ + im^/h imj/h, (23) 

q « {zm7/?i}i/% (24) 

C_, i?2 0, (25) 

C+^l. (26) 

In this case the state vector for a; > becomes simply the plane wave with wave number k in the ground state. This 
means that for increasing 7 there is less and less reflection, but also less and less absorption, i.e. photon detection, so 
that the laser has less and less effect on the atom. 

At first sight the occurrence of reflections may seem surprising since the laser only couples to the internal degrees 
of the atom and since only applies to the time development before the first photon detection. Physically this can 
be understood from the coupling of the atom to the quantized electromagnetic field. The laser changes the internal 
state, this in turn changes the quantized electromagnetic field and its momentum distribution. This in turn changes 
the momentum distribution of the atomic motion. Mathematically the reason is of course the step function in front of 
the matrix, similar as for a square-well potential. The consequences of the nonzero reflection will be discussed further 
below. 

By decomposing an initial state as a superposition of eigenfunctions one obtains its conditional time development. 
This is easy for an initial ground-state wave packet ('^(^'*'')) coming in from the far left side and with to in the remote 

past. Indeed, if ipik) denotes the momentum amplitude the wave packet would have as a freely moving packet at 
t = 0, then 

*(x,i)=/ dkiP{k)^k{x)e-'^^'''^"' (27) 



describes the conditional time development of a state which in the remote past behaves like a wave packet in the 
ground-state coming in from the left. 
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III. THE REFLECTION PROBLEM AND THE NO-DETECTION PROBABILITY 

From Eq. (||) one obtains, with ^ — Clj^J)), 

N, = y'dx{|V'W(a;,i)|2 + |^(2)(x,t)|2} (28) 

and Eq. ph becomes 



U{t)^^{m\Hc-HMt))- (29) 
Since He — H] = —i^h\2){2\^ the first-photon probabihty density is given by 

/oo 
dx\ilj^''\x,t)\'' . (30) 
-OO 

The probabihty of no photon detection at all is, for to in the remote past, 

/oo 
dt'n(i') = A^t=oo. (31) 
-oo 

For physical reasons, only i/;*^^) contributes to this, and only for a; < 0. The latter follows from the fact that, for a: > 0, 
the ground-state part will eventually be pumped by the laser to the excited state. Since t = oo in Eq. (|3l|), only the 
reflected part remains, and hence the no-detection probability becomes 

/>oo 

Wt=oo= / dk\R^{k)\^\i^{k)\^. (32) 

As a consequence, n(t) is in general not normalized to 1. Physically the probability for missing an atom increases 
with VL, the strength of the laser driving. This is also seen mathematically from the expression for i?i in Eq. (|T^). 
An example is given in Section ^ (cf. Fig. ^ further below) and a practical approach to bypass this problem is also 
discussed there. 

On the other hand, for fc — > cx) reflection becomes negligible since then Ri 0. Hence for faster atoms reflection 
does not pose a problem. This will also be exploited for practical purposes in Section 

IV. DELAYS VERSUS REFLECTIONS AND AN IDEALIZED DISTRIBUTION 

The approach to quantum arrival times by means of first-photon detections contains a built-in "delay" since an 
excited atom will not emit a photon immediately, due to the finite decay rate 7. This is in addition to the time the 
photon takes to reach the detector; the photon travel time, however, can easily be taken into account, so it will not 
be considered here any further. 

It seems natural to try to obtain an ideal arrival-time distribution by considering faster and faster decay times, i.e. 
taking 7 — > 00, at least theoretically. This will not work, however. The reason is that for increasing 7, with all other 
parameters kept fixed, the driving by the laser becomes less efficient so that, in the limit 7 — > c», the wave packet 
remains unaffected, with no excitation and no reflection, as can be seen from Eqs. (|2^) - ( p6|) above. Moreover, if both 

and f2 go to infinity with 7/fi kept fixed, then Ri —t —1 and everything is completely reflected without excitation 



One might also be tempted to avoid reflection by choosing weak driving, $1/7 ^ 1. This, however, would cause 
a severe delay problem since the laser would take more time to pump the atom to the excited state. Hence the 
first-photon emission would also take more time. To see how relevant detection delays are, we have compared H(t) 
with the flux and the axiomatic probability distribution HK(i) of Kijowski for a Gaussian wave packet. The result 
is given in Fig. |l|. Depending on the parameters, the delay and reflection problem may be either very relevant or 
negligible. A detailed analytic investigation of this question is given in Section 0. 

A way out of the conflicting problems of reflection (missed atom) and increasing delay times for weaker driving is 
the transition from the "experimental" H(t) to an idealized arrival-time distribution, obtained as follows. A two-level 
atom at rest, when driven by a resonant laser, has a deflnite probability density, W{t), for the detection of the flrst 
photon, given by ||2^ 

W{t) = ^e-^*/2|e^*/2 - e-^*/2pe(i), (33) 
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FIG. 1: Time-of-arrival distributions: Flux J (solid line, here indistinguishable from Kijowski's IIk) and D (first photon, dots). 
Note the delay in H. The initial state is a minimum-uncertainty-product Gaussian for the center-of-mass motion of a single 
Cesium atom in the ground state with (v) — 9.0297 cm/s, (x) = —1.85 fj,m, and Aa; = 0.26 fj,m.; Q = 0.09997; s-U figures are 
for the transition 6^P3/2 - &^Si/2 of Cesium with 7 — 33.3 MHz. 



with 

5= 1(7^ -4f72)i/2. (34) 

Intuitively, the delay-time mechanism for a moving atom ought to be similar to that for an atom at rest. Should it 
then not be possible to somehow compensate the delay in Il{t) by that of the atom at rest and thus arrive, in some 
limit, at a delay-free ideal distribution? To achieve this, we assume the (experimental) arrival-time distribution Tl{t) 
to be the convolution of a hypothetical ideal distribution, Hid, with the distribution W{t) for an atom at rest, 

n = Hid * W. (35) 

The delay in 11 is then mainly due to that contained in W. 

The hypothetical ideal distribution Ilid is obtained by a deconvolution via Fourier transform from Hid ~ H/W 
where U{v) = J dte-'^^Uit) etc. From Eqs. {M) and iM) one finds 



U{iy) =-fJdxJdkJ dk'^{k)i}{k')(j)^^\x)c^^,\x)2TrS - ^{k^ - fc")) (36) 

and from Eq. (||) 

Wi.,) ^ ^ 5!2Zi , (37, 

and thus 

I + Cxiu + C2{wf + C:i{iuf , (38) 



W{v) 
where 
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In the time domain this gives 

nid(t) = n{t) + Ciii'it) + C2ii"{t) + C3n"\t) . (40) 

In the hmit of no reflection, for 7 —>■ 00, the delay problem pointed out above for 11 should be absent for Hid- In 
fact, in this limit one sees from Eqs. ( pc| ) - ( p^ ) and from Eq. (|3|) that in l/W only iv^jV? remains relevant for 
large 7. Furthermore, i?2 ^ for 7 ^ cx), and since e"*''^ — > in Eq. (^6|) the integral over a; < goes to zero. 
Similarly the terms in (\)''^^ (x) and <^^y (x) containing e~**^-^ and g''^'-^ drop out. For 7 ^ 00 one has, from Eqs. ( [20| ) 
and (p^), C+ ^ 1 and the integration of exp{— ifc+x + i/c^x} over x > gives i/{k'^ — Hence 



By the (5 function, v = h{k — k')(k + k')/2m. Hence in the limit 7 ^ 00 one obtains, by Eqs. ( |20| ) - (p^), 

£^ [ dkdk'J{k)M')(k + /c')27r(5 fi/ - — (fc2 - fc'2)^ . (42) 

W{v) 2TT2m J V 2m y 

Going to the time domain one finally obtains 



1 h 



271 m I 2 



h 

2mi 



T{^(0,i)V'(0,i)-^'(0,i)V(0,t)}, 



which is the flux for the free wavefunction 'ip{x,t) at a; = 0, i.e. without laser [ p2[ . 

This is an extremely interesting result since is a natural candidate for the arrival-time distribution, as pointed out 
in the Introduction. We note that J,^, is normalized to 1 for a particle which has only positive momentum components. 
This is seen for example from Eq. (^) for v — 0. 

The limit in Eq. (^3|) means that Hid can be approximated by for sufficiently large 7. Physically, it is important 
to determine the parameter ranges for which this approximation is a good one. For this to be valid a simple sufficient 
condition on the parameters can be derived as follows. First one considers the case ri^/7^ <C 1 (weak driving) and 
makes the corresponding approximations in Hid(i). Denoting by AE a measure of the magnitude of 

|l|fc2-fc'2|, (44) 

as determined by V'('s) and ijj{k'), one then considers the case 7 3> AE/h and AE/h ^ i^^/7- Then one obtains that 
Hid is close to J^. The inequalities can be written in the form 

n'^/j < AE/n < 7 (45) 

Thus Hid can be replaced by if these inequalities are satisfied. Even outside this parameter range, may be an 
excellent approximation to Hid, as shown in Fig. ^. The convergence of Hid to shows that also Hid may contain 
small negative values. Since this occurs neither for H nor for W, the intuitive ansatz of Eq. ( |35| ) cannot always be 
fulfilled with a strictly positive distribution. The reason for this clearly is that the ansatz of a convolution in Eq. ( |35| ) 
is too simple and ought to be replaced by something more sophisticated. On the other other hand, this result gives a 
handle at the quantum mechanical probability flux and indicates a method how to measure it. 

The above deconvolution procedure which recovers the quantum mechanical particle flux essentially works because 
the weak-excitation limit taken allows a clear separation of (1) the time dependence associated with the motion of the 
wavepacket, and (2) the time dependence associated with the internal degrees of freedom (excitation, Rabi oscillation, 
and decay). It seems reasonable that this might be done. However, the weak-excitation limit implies that the waiting 
times for the first scattered photon are of the order j/fl'^ and therefore very long. Hence the number to be measured 
is essentially to be obtained from the subtraction of two very large numbers. Experimentally this is a difficult thing 
to do with high accuracy and requires small measurement errors. For practical purposes it is therefore important to 
know when delays and reflections can be safely neglected, since then the transition to Hid by deconvolution is not 
necessary. This is investigated quantitatively in the next section. 



8 




FIG. 2: Excellent agreement between Hid (filled circles) and J (solid line); deviations from Hk (dotted line) and 11 (dot-dashed 
line). The initial wave packet is a coherent combination ip = 2~^^'^(il>i + 1P2) of two Gaussian states for the center-of-mass 
motion of a single Gesium atom that become separately minimal uncertainty packets (with Axi = Ax2 = 0.021 fim, and average 
velocities {v)i = 18.96 cm/s, {17)2 = 5.42 cm/s) at s = and t — 2 /is; SI = O.377. 



PARAMETER RANGES WITH NEGLIGIBLE DELAY AND REFLECTION. 



In the case of weak driving, il <C 7, the excited state population is neghgible compared to that of the ground state. 
Hence the reflection coefficient for the ground state is the only one that matters in this case, 



\R,{k)\^^^{n/jrm/Er, 



(46) 



This is small ii E ^ hfl. For strong driving, fl ^ j, the reflection coefficients take a simple form when the energy E 
of the plane wave satisfies 

E > nn. (47) 
Then both states are populated roughly equally for a; > 0, and with Eq. (^) the reflection coefficients become 

(48) 



Both coefficients are small if Eq. ( ^ ) holds. 

To quantify the detection delay, let us define as the difference between the average time of the first photon 
emission, (i)n = dttll{t), and {t)j = dtt we note that the "average arrival time" at a; = evaluated 

with the flux at a; = 0, {t)j coincides with the average of Kijowski's distribution ||l|, For negligible reflection, 
'4'(a:,<) is, for a; < 0, nearly the same as the free wavefunction, and therefore, by a partial integration. 



/oo 
dt (for negligible reflection) , 

-00 



(49) 



where 7V+ = /(,°°(|i/;(i) (a;, i)|2 + |,/>(2) (x, t)|2) dx. Hence, in the case of negligible reflection the delay is associated with 
the amount of penetration of the wave into the laser region. For weak and strong driving a straightforward calculation 
yields the simple results 



Td « 7/^^^ (weak driving) . 
Td w 2/7 (strong driving). 



(50) 
(51) 
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FIG. 3: Negligible delay and reflection with Eq. (|53|): 11 (first photon, dots) and (solid line, here indistinguishable from 
Hk) for initial Gaussian parameters (v) = 90.30 cm/s, (x) = —218.02 ^m. Ax — 26.46 ^m; f2 = 57. 

It is interesting to note, and physically very reasonable, that this coincides with the average time between two photon 
emissions for a two-level atom at rest, driven by a resonant laser, as easily seen from Eq. (|33|). If At denotes the 
width of n(i), then one needs ^ At for the delay to be negligible. 

If one denotes by E the average energy for wave packets that are sharply peaked in energy or, more generally, as 
the smallest significant energy component, then the above conditions for negligible reflection and negligible delay can 
be summarized as 



Fig. y shows a striking example for which these conditions are fulfilled. Here the first-photon distribution 11 is 
indistinguishable from and from Kijowski's distribution. 

For strong driving the delay due to the time needed for pumping the atom to an excited level is small, but on the 
other hand there is a large probability to miss out the atom altogether if the condition E ^ hfl is not fulfilled. As a 
consequence, the first-photon probability density n(t) will not be normalized to 1. Let us then define a normalized 
distribution Il]\j(t) — Tl{t)/ J dt' Il{t'). For practical purposes this can be amazingly efficient in some parameter 
domains as Fig. |j shows. There, 11 is far from being normalized to 1, but IIat and J^, which in this example has no 
negative parts, coincide beautifully. Interestingly, one notices a small difference to Kijowski's distribution. It might 
also be possible to use the general normalization procedure in terms of operators proposed in Ref. 



In this paper, we have investigated a proposal to determine arrival times of quantum mechanical particles. The 
proposal is based on the intuitive idea to illuminate the arrival region by a laser and to consider a traveling two-level 
atom. The time of the first emitted photon is then taken as a measure for the arrival time. By repeating the experiment 
one obtains a probability density, n(t), for the time of the first photon. We have discussed for the one-dimensional 
case in what way Il{t) can be regarded as an atomic arrival-time distribution. Restrictions arise from refiections and 
delays. Reflections originate from the interaction with the laser and delays from the time needed for the pumping 
and the ensuing photon emission. The natural idea that an ideal or an axiomatically proposed distribution might 
be obtained from Il{t) in the limit of a very strong or very weak laser and very large Einstein decay coefficient of 
the excited level has turned out not to be true. However, and this is a main theoretical result of the paper, one 
can subtract the delay by a deconvolution with the first-photon probability density for an atom at rest and then, 
surprisingly, for larger and larger Einstein coefficient one obtains the quantum mechanical probability current J as 
the limit distribution. This quantity J has previously been considered on axiomatic grounds as a candidate for the 
arrival time distribution, and the connection of H(i) with J indicates, to our knowledge for the first time, a way 



Ti/E < -ijVp- < At (weak driving) 

YijE <iCVr^ < 2/7 < At (strong driving). 



(52) 
(53) 



VI. CONCLUSIONS 
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FIG. 4: Improvement by normalization: IIjv (white circles) right on top of (dashed line); also shown 11 a' (solid line) and 
n (filled circles) . The initial Gaussian wave packet is chosen to become minimal when its center arrives at x = (in the absence 
of the laser) to enhance the difference between Hk and J^; (v) = 0.9 cm/s, Aa; = 0.106 fj,m, Q = 87. 



to measure the quantum mechanical probability flux. We have also determined parameter domains for which the 
deconvoluted expression is already sufhciently close to J. Although the non-deconvoluted Il{t) is not the same as J 
and the axiomatically proposed distribution of Kijovirski, it can, for experimental purposes, approach the latter two 
sufficiently closely. Parameter domains for which this holds have been explicitly determined; this is another main 
result of the paper, more of a practical nature. 
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